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Abstract 

We use finite-size scaling of Lee- Yang partition function zeroes to study the critical be- 
haviour of the two dimensional step or sgn 0(2) model. We present evidence that, like the 
closely related XY-mode\, this has a phase transition from a disordered high temperature phase 
to a low temperature massless phase where the model remains critical. The critical parameters 
(including logarithmic corrections) are compatible with those of the XY-model indicating that 
both models belong to the same universality class. 
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1 Introduction 



In a recent paper [|J, we demonstrated the power of finite-size scaling applied to Lee- Yang zeroes 
[g] in uncovering logarithmic corrections to scaling in the two-dimensional XY- or 0(2) spin model. 
In this paper we apply the same techniques to the closely related 'step model' [||, also known J|, ||] 
as the sgn 0(2) model. The question of criticality of this model has, until now, been unresolved 
despite several analyses based on high temperature series (see || for a review) and on numerical 
simulation ffl, B. The interest in the model arises from its possible membership of the XY-model 
universality class which exhibits the Kosterlitz-Thouless [KT] phase transition ||. Like the XY- 
model, the step model has a configuration space which is globally and continuously symmetric. 
Unlike the XY-model, however, the interaction function is discontinuous and the Mermin-Wagner 



theorem [ 10 1 does not apply. Nonetheless, it is expected that if a phase transition exists in the step 



model, it should not be to a phase with long range order [|], ||, |TTJ. 

Sanchez-Velasco and Wills || presented evidence of critical behaviour starting at (3 C = 1/T C = 
0.91 ± 0.04. This was based on finite-size scaling [FSS] of the spin susceptibility. Since the as- 
sociated critical index r){T c ) was significantly greater than that measured for the XY-model, it 
was concluded that the step and XY models are not in the same universality class. In this paper 
we present evidence that the step model is not critical at that temperature. However it is critical 
at lower temperatures with a critical index rj(T c ) compatible with the XY value. The accuracy 
afforded by the Lee- Yang zeroes study is a crucial part of the analysis. 



2 The step model and the XY— model 

Consider the partition function 

Z(j3,h) = Y,e~ f>H+Mi '", (1) 

where the Boltzmann factor is (3 = 1/kT, n is a unit vector defining the direction of the external 
magnetic field and h is a scalar parameter representing its strength. The summation is over all 
configurations open to the system and s x is a unit length two-component spin at each site x in the 
cubic lattice A = L d (d = 2). The magnetisation for a given configuration is 

M = 4 • 

In the case of the XY-model, the interaction hamiltonian is 

d 

HxY = — X! X! ^ x ' Sx+n , 

xeA/i=i 

while for the step model it is 

d 

Thus the continuous (cosine) dependence of the interaction energy in the usual XY-model is 
replaced by a discrete step function dependence. 

The leading infinite volume critical behaviour of the 2D XY-model is characterized by exponential 
divergences (essential singularities) in the thermodynamic functions ||. In terms of the reduced 



temperature t = 1 — (3/ ' [3 C — » + the (leading) infinite volume scaling behaviour of the correlation 
length and the zero-field magnetic susceptibility is (respectively) || 

&o(t) ~ e ai "\ (2) 

xoo(t) ~ e^, (3) 

where = 1/2 and r/ = 1/4. The .YY-model remains critical for all (3 > (3 C . 

For models obeying the Lee- Yang theorem @, the partition function zeroes in the magnetic field 
strength (h) plane (the Lee- Yang zeroes) are all on the imaginary axis. In the high temperature 
phase these zeroes remain away from the real axis, pinching it only as t — > + (in the thermodynamic 
limit). The zero closest to the real axis marks the edge of the distribution of zeroes and is known 
as the Yang-Lee edge, zyl- The theorem has been proved only for certain models, the XY-model 
included [l^] . In [0 we used this fact to show the above leading critical behaviour for the X Y-model 
and the corresponding behaviour of the Yang-Lee edge, zyl are, in fact, modified by logarithmic 
corrections: 



>oo L 

z Yh(t) ~ (5) 



Xoo(t) ~ & v t\ (4) 

-A 



where 

1 1 5 

X = --(d + 2- V ) = -^ (6) 

and the parameters r and p (= — r/2) are logarithmic correction indices Q. The corresponding 
FSS behaviour for the susceptibility and the first zero z\ (= zyl) at t = is jlj 

Xl(0) - L 2 ^(lnL)-^, (7) 
2i (0) ~ L A (ln L) r . (8) 

For the 2D XY-model the numerical value of r was found to be small but non-zero (— 0.023±0.010) 

The objects of the present analaysis were to establish 

1. if the scaling behaviour of the very precisely determined Yang-Lee edge would unequivocably 
determine whether the step model had a critical phase 

2. if so, whether the phase transition is in the same universality class as the XY-model. 



3 Method and results 



The methods used are those described in |lj. A single cluster algorithm [ 13 1 is used to generate a 
large number of measurements (100K for each lattice size L and temperature 1//3) of the energy 
H and the magnetisation M at zero external magnetic field. Such data were obtained for lattice 
sizes L = 32, 48, 64, 128 and 256 covering the (3 range 0.86 to 1.40 with varying degrees of spacing. 
Multi-histogram techniques [jlj, [H| were used to combine data at different values of j3 and so 
obtain detailed (3 dependence. In the neighborhood of possible critical points we used sufficiently 
fine spacing (typically 0.025) to ensure adequate overlap between histograms for a given size of 
lattice. 

The partition function for a complex magnetic field (h = h r + ihi, h r ,hi € K) can be written in 
terms of a real and an imaginary part 0, [l6|, [l7| as 



Z(P, h r + ihi) = ReZ(/3, h r + ihi) + iImZ(/3, h r + ihi) > 



(9) 



where 



ReZ((3,h r + ihi) = Z (P,h r ) {cos (hiM)) p>hr , (10) 
hnZ(p,hr + ihi) = Z (13, h r ) {sin (hiM)) p , hr . (11) 

Here the subscripts indicate that the expectation values are taken at (inverse) temperature /3 and 
in a (real) external field h r . 



No specific proof of the Lee- Yang theorem exists for the step model. We can, however, use (|10|) and 
( |IT|) to (numerically) determine the loci along which the real and imaginary parts of the partition 
function separately vanish. These formulae concern expectation values of real quantities in a real 
external field and at no stage is a simulation involving a complex action involved. The Lee- Yang 
zeroes are then the points in the complex /i-plane where the loci intersect [17]. We have determined 



these loci and thereby the Lee- Yang zeroes close to the real /i-axis. We were able to determine the 
first 15 zeroes and found that they lie on the imaginary /i-axis for all the lattices studied. Thus we 
have numerical evidence that the step model obeys the Lee- Yang theorem. The following analysis 
applies to the first zeroes only (the Yang-Lee edge) and we defer the study of higher zeroes to a 



later paper [18| 



The analysis began with a rough search for the leading critical behaviour predicted by (||). An 
independent test was also made using the (less accurate) susceptibility data and (0). Both methods 
indicated critical behaviour setting in for (3 ^ 1.2. In Fig. |l| we show a typical log-log plot of 
z\ the Yang-Lee edge versus L. This is at a typical candidate value of the critical temperature 
((3 = 1.22). The errors are considerably smaller than the symbols. For example, at (3 = 1.22 we 
found zx = 0.0024136(7) and 0.00017902(6) at L = 32 and 128 respectively. The slope of Fig. | 
gives the effective leading index A g g- ignoring corrections. At (3 = 1.22 this is X q q = —1.8761(2) 
where the chi-squared per degree of freedom (x 2 /dof) for the linear fit shown is 0.85. In Fig. 2(a) 
we display the result of such fits as a function of (3. The effective exponent A g g- is just the slope 
of the log-log linear fit which should obtain if critical behaviour is present (ignoring logarithmic 
corrections). The corresponding x 2 /dof for the linear fit is also shown. Acceptable values are only 
found for j3 in excess of around 1.2. To quantify this statement we demand 

X 2 /dof < 2.0 (12) 

which means (3 > 1.185. We note that the corresponding values of A e fj ( < — 1.872(2)) include 
that (—15/8 = —1.875) corresponding to the KT prediction. 

Fig. § is evidence for (i) the validity of FSS over a range of values (3 > (3 C ~ 1.185 and (ii) at j3 c the 
exponent A g g very close to the expected KT value (-15/8). Observation (i) means that, as in the 
XY case, the system remains critical for (3 > (3 C and (ii) is evidence that A is in fact -15/8 and the 
question of whether or not the step model belongs to the same universality class as the XY model 
must now be answered by determination of the correction exponent r. 

We therefore assume the behaviour (|8]) with A = — 15/8 at = (3 C . The expected leading behaviour 
is removed and linear fits to 

ln(ziL 15/8 ) vs. InlnL (13) 

performed. The results are shown in Fig. 3. Since the value A = — 15/8 (r] = 1/4) is only expected 
at (3 C these results can be used to identify the possible values of critical temperature and to test for 
the presence of logarithmic corrections as in the X Y-model ||] . 



Applying the same criterion ( |12|) as for the leading behaviour, we search for a range of /3 C values 
giving an acceptable fit. We find 

1.195 < (3 C < 1.295 and correspondingly, 0.009 > r > -0.034 . (14) 



The range of acceptable r values includes that found Q for the YY-model (—0.023 ± 0.010) and 
that corresponding to no logarithmic corrections (r = 0). As in our previous work Q, this range 
excludes the prediction r = —1/16 coming from an approximate renormalisation group treatment 
of the XY-model |J. Thus we conclude that the present data are compatible with the step 
model being in the same universality class as the YY-model. We do not, however, exclude other 
possibilities. 

The susceptibility data are consistent with the above analysis. If one assumes the KT value of 



v(Pc) = 1/4, one can construct a so-called Roomany-Wyld beta function approximant [19] from 
the finite-size data and use its zero to locate (3 C These approximants, based on pairs of lattice 
size L,L', converge very rapidly |L9|]. We estimate j3 c = 1.22 ± 0.02. This last analysis of course 
neglects possible logarithmic corrections to scaling. 

We have also studied the specific heat. As for the XY-model, the step model data show a broad 
peak with no obvious relationship to the position of the leading critical point. The finite-size 
dependence is not dramatic and is likely to be of little value in further elucidating the criticial 
behaviour. A related question is to what extent one can make use of the Fisher zeroes 
i.e. zeroes of the partition function in the complex (3 plane at zero external magnetic field h. For 
both this and the YY-model, these are much harder to locate than the Lee- Yang zeroes and are 
consequently less accurately determined. 



4 Conclusions 

The use of finite-size scaling applied to Lee- Yang zeroes has allowed us to present detailed evidence 
of critical behaviour in the two-dimensional step (sgn 0(2) spin) model. The data are consistent 
with this model being in the same universality class as the XY- ( 0(2) spin) model. That is, 
it undergoes a Kosterlitz-Thouless type transition with susceptibility index r]([3 c ) = 1/4 and we 
determine that (3 C lies in the range 1.195 < (3 C < 1.295. With the available statistics, we found the 
logarithmic correction exponent to lie in the range —0.034 < r < 0.009. This should be compared 
with our measurement for the XY correction exponent M, —0.033 < r < —0.013 with which it is 
compatable. The step model results are however also compatable with no logarithmic corrections 
(r = corresponds to j3 c = 1.21 and x 2 /dof = 0.92). 

The Mermin- Wagner theorem [1C] does not apply directly to the step model because of the dis- 
continuous nature of the interaction hamiltonian. However, it has long been believed [^, H that 
if a phase transition does exist, it will not involve a phase with long range order. The evidence 
presented here supports this view. 

This raises a question as to the nature of the mechanism driving the phase transition in the 
step model. The KT phase transition of the YY-model is understood to be driven by the bind- 
ing/unbinding of topological solutions (vortices). However, the energetics of vortex formation are 
very different in the step model [11, |f| H|. Since vortices with effectively zero excitation energy can 



be created at all non-zero temperatures, the usual KT argument does not naturally lead one to 
expect such a phase transition in the step model. 

If this is indeed the case, some other driving mechanism must be responsible for any phase transition. 
It would then be remarkable if — as the evidence presented here indicates - both models belong 
to the one universality class. 
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Figure 2: Test of expected leading critical behaviour: (a) the effective exponent X c g (slope of a 
straight line fit) and (b) the corresponding x 2 /dof versus (3. 
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Figure 3: Logarithmic corrections: (a) the logarithmic correction exponent r to the Yang-Lee edge 
is shown as a function of the assumed critical coupling [3 C and (b) the corresponding % 2 /dof. 



